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Abstract 

The structures of compressed rubidium polyhydrides, RbH n with n > 1, and their evolution 
under pressure are studied using density functional theory calculations. These phases, which start 
to stabilize at only P = 2 GPa, consist of Rb + cations and one or more of the following species: 
H~ anions, H2 molecules, and units. The latter motif, the simplest example of a three-center 
four-electron bond, is found in the most stable structures, RbHs and RbH3. Pressure induces 
the symmetrization of H^~. The thermodynamically most stable polyhydrides metallize above 
200 GPa. At the highest pressures studied, our evolutionary searches find an RbH6 phase which 
contains polymeric (H^ - )^ chains that show signs of one-dimensional liquid-like behavior. 
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I. INTRODUCTION 



H3" has been known to mass spectroscopists since Thomson reported it in 1912 [lj, but 
it took 90 more years to reliably observe H3 [2HI]. Already in 1935 Coulson predicted the 
triangular configuration of the cation, and the linear arrangement of its anionic cousin [5]. 
Theoretical studies confirmed these geometries [6] but also concluded that H3 was unstable 
with respect to decomposition into H 2 and H~ [7J. In the 1990's ab-initio computations 
became reliable enough to show that H3 is a minimum on the potential energy surface 
(PES) [HlEj, which has been subsequently explored [101 HI]. We came across the H3" anion 
serendipidously during our computational quest to predict hydrogen-rich structures that are 
metallic and potentially superconducting at experimentally accessible pressures pT2l [T3] . 

Hydrogen turns black at 320 GPa [14], but it does not become metallic even at pres- 
sures approaching those found in the Earth's core [15j. Once metallized, H 2 is predicted to 
become superconducting at temperatures as high as 240 K [16]. On the other hand, band 
gap closure and concommitant superconductivity in various metal hydrides such as GeH 4 
[T7] , SnH 4 [18] , AIH3 [19l [20] , and GaH 3 [21] has been computed for pressures attainable 
in diamond anvil cells. Various attempts to dope hydrogen by an impurity cab able of pro- 
moting metallization have been carried out [22] . Compression not only affects the electronic 
structure and properties of extended systems, but it also modifies their chemistry: it may 
entice elements which would not normally combine to do so [23] , or to mix in non-traditional 
proportions. An example of the latter is the van der Waals compound Xe(H 2 )7 [24], or the 
SiH 4 /H 2 phases with ratios of 1:1, 1:2, and 1:5 j25j [26] which have been observed in recent 
experiments. 

The exciting findings in high-pressure research have inspired us to theoretically explore 
the potential energy landscape of compressed alkali metal poly hydrides, MH n with n > 
1. These stoichiometeries are unstable with respect to decomposition into MH and H 2 at 
atmospheric conditions, but calculations show that the enthalpies of formation of the most 
stable lithium [27] and sodium polyhydrides [28] become negative at about 100 and 25 GPa. 
The pressures necessary to stabilize these systems were found to correlate with the ionization 
potentials (IP) of the metals, and the preferred stoichiometries (LiH 6 and NaH 9 ) with the 
size of the alkali metal cation. 

Rb has an even lower IP than Na, and Rb + a larger radius than Na + . In order to 
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determine how this affects the rubidium poly hydrides, we have performed first-principles 
calculations coupled with the evolutionary algorithm (EA) XtalOpt [29j [30] to search for 
the most stable RbH n up to 250 GPa. The results of our studies on RbH n (n > 1) illustrate 
that under pressure: (i) compounds with unusual stoichiometries may be stabilized; (ii) 
there is an increased tendancy to form multicenter bonds; and (iii) orbitals which chemists 
do not usually think of (ie. Rb 4<s and 4p) become important. 

Below we explore the various phases of the rubidium polyhydrides at pressures ranging 
from 2 GPa — the onset of stabilization — up to the eventual insulator-to-metal transition 
(MIT) which occurs just over 200 GPa. RbH 5 , the most stable structure throughout a large 
part of the pressure range considered, contains H 2 molecules, Rb + cations and H3 anions 
which become symmetric above 30 GPa. Other fascinating phases, such as RbH 9 and RbH 3 
also have domains of thermodynamic stability. At the greatest compressions considered, the 
hydrogen sublattice itself shows signs of polymerization. 

II. RESULTS AND DISCUSSION 

The calculated enthalpies of formation, AH F) of the most favorable RbH n structures 
found in our evolutionary runs are provided in Fig.[j](a). RbH 2 (33.3 mol% H 2 ) is not shown 
since it was the only phase found to be unstable with respect to decomposition into RbH 
and H 2 at every pressure considered. At only ^2 GPa RbH n with n=5-9 become preferred. 
This industrially attainable pressure is an order of magnitude lower than that necessary to 
stabilize the sodium polyhydrides, highlighting again the import of the alkali metal's IP for 
the formation of these phases. 

If a tie-line is drawn connecting AH F (A) with AH F (B) and AH F (C) falls beneath it, 
then A and B will react to form C. If another phase D falls above the tie-line connecting A 
and B, it is expected to decompose into these two structures, provided the kinetic barrier 
is not too high. The tie-lines may be used to construct a convex hull — a point lying on 
it corresponds to a thermodynamically stable phase at the given pressure [23]. The dashed 
line in Fig. [j](b) traces out the estimated convex hull of RbH n at 100 GPa, where RbH 5 
has the most negative enthalpy of formation. Other thermodynamically stable phases are 
RbH 3 , RbH 9 , and RbHi 6 . Above this pressure the AH F of all of the phases — with the 
exception of RbH 3 — increases. The destabilization of structures at greater compression was 
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FIG. 1: (a) The enthalpy of formation, AHp, for the reaction RbH+ ^(H-2)n-i RbH n versus the 
H2 composition at varying pressures, where the enthalpy of RbH and H2 are computed for the most 
stable structures from Refs. [31] and [32]. The corresponding convex hulls at (b) 100 GPa and (c) 
250 GPa are shown below to highlight thermodynamically stable RbH n structures. The horizontal 
axes in (b) and (c) are not to scale. Filled symbols represent points at which EA structure searches 
were explicitly run; open symbols represent structures which were re-optimized at said pressures 
from an EA search performed at another pressure. 



not previously observed in studies of the lighter alkali polyhydrides up to 300 GPa [271 EE] • 
At 250 GPa RbH 3 is the most stable structure, with RbH 6 , RbH 9 and RbHi 4 falling on the 
convex hull. Comparison of these results with previous findings for LiH n and NaH n suggests 
that there is no simple relationship between the preferred stoichiometries of the alkali metal 
polyhydrides and the radius of the cation. 

Within the pressure range considered, the most stable phases contain between 3 and 9 
hydrogen atoms per rubidium. They may all be reasonably divided into four groups which 
are based on qualitative assessments of the atoms' spatial arrangements: 1) Rb + cations 
and (H 2 ) -n / 2 molecules with bonds longer than 0.75 A; 2) Rb + , H 2 molecules, and hydridic 
H~ anions; 3) Rb + , H 2 and H3 molecules; and 4) Rb + and H3 . The first two structural 
makeups have already been observed in LiH n and NaH n . The latter two groups have not 
yet been noted, but are found in the most stable rubidium polyhydrides, RbH 5 and RbH 3 . 
RbH 6 consists of arrangements of H3 units brought together to form what are essentially 
chains. Herein, we examine the geometries of these phases, and their electronic structures 
— especially the imminent metallization under compression. 

A. Pressure— Induced Stabilization of 

Below 100 GPa the most stable stoichiometries are RbH 9 and RbH 5 . The enthalpy of 
formation of a Pm-symmetry RbH 9 phase becomes negative near ^2 GPa and RbH 9 remains 
the preferred combination up to ~16 GPa, at which point it is overtaken by RbH 5 as shown 
in Fig. [2](c). The lowest enthalpy structures differ at 10 and 100 GPa, so below we refer to 
them as the "low pressure" (lp) and "high pressure" (hp) structures. 

The linear motifs which characterize the H3 units within the Cmcm-symmetry hp-RbH 5 
phase are highlighted in Fig.[2](a) [54J. At 100 GPa, the intramolecular H-H distances in the 
H3 and H 2 units measure 0.91 A and 0.77 A. The latter is slightly elongated when compared 
with the 0.73 A in pure P63/m-H 2 j32] at this pressure. The shortest intermolecular H- 
H distances in hp-RbH 5 are roughly equivalent to those in pure H 2 (1.56 A vs. 1.50 A) 
implying that RbH 5 can reasonably be thought of as containing an arrangement of H 2 and 
H3 molecules. 

From visualization of the phonon modes of hp-RbHs lying above 2000 cm -1 , they are 
assigned to H3" molecular vibrations between 2000-2700 cm -1 and H 2 near 4000 cm -1 . 



5 



These modes are shown schematically in the phonon DOS in Fig. |2j(d), and could serve as 
a fingerprint to characterize RbH 5 , particularly the symmetric H3 stretch near 2500 cm -1 . 

The symmetric H3 units in hp-RbH 5 persist as the pressure is decreased until 15 GPa, 
wherein the system becomes thermodynamically unstable with respect to lp-RbH 5 with 
asymmetric H~ • • • H 2 units. Actually, at around this same pressure both RbH 5 phases are 
predicted to decompose into RbH and RbH 9 . The lp-RbH 9 shown in Fig. [2](b) at 10 GPa 
also adopts asymmetric H~ • • • H 2 arrangements but, unlike lp-RbH 5 , each H~ atom has 
several weak (linear) H~ • • • H 2 interactions. At 10 GPa, the four closest H~ • • • H 2 contacts in 
lp-RbHg measure ^1.85 A but there is only one such contact in lp-RbH 5 measuring 1.27 A. 
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FIG. 2: 2x2x2 supercells of the (a) hp-RbH 5 phase at 100 GPa, and the (b) lp-RbHg phase 
at 10 GPa. The hydrogen atoms with H~ or character are shaded to enhance visibility, (c) 
AH F = H(RbR n ) - (^)if (H 2 ) - #(RbH) for these structures, (d) Phonon calculations of RbH 5 
at 30 and 100 GPa did not reveal any imaginary frequencies, indicating that these are local minima. 
The phonon DOS for hp-RbH 5 at 100 GPa is shown. 
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The H~ • • • H 2 arrangements adopted in these structures conform to the expectation of 
an increase in multicenter bonding in compressed solids [33] . Such a precedence has already 
been set in pure H 2 wherein calculations predict that at a few TPa triangluar H^j" motifs may 
form along with a background of negative charge [34] . The electronic structure of H^j" and 
H3 can be understood in terms of three-center bonding schemes; for two electrons (3c-2e) 
electron poor, and four electrons (3c-4e) electron rich. 

A classic example of a 3c-4e bond is the triiodide anion. In fact Rb forms one of the 
salts of I3 [35J which, if you consider Rb + H3 as an alkali trihalide, also agrees with the 
observation that pressure can favor configurations adopted by heavier elements in the same 
group [33] . 

The shallow energy minima in the H3" PES yields H-H distances (d and z depicted in Fig. 
o measuring: 0.75 and 2.84 A [TO]. The lowest energy barrier connecting the two minima is 
a symmetric H3 molecule with H-H distances of ^1.06 A [10]. Such symmetric H3" motifs 
have also been discussed in the literature as transition states in hydrogen exchange processes 
of M(H)(H 2 ) metal complexes [36] . 

The observation of symmetric H3" units at 100 GPa and their asymmetrization as the 
pressure is lowered (~10 GPa) suggests that the equalization of d and z is pressure-driven. 
To test this further a computational experiment was carried out wherein the total length 
(d + z) of an H~ • • • H 2 unit was constrained to a fixed value in order to mimic an applied 
external pressure, and the molecule was subsequently optimized [55] . In Fig. [3] we give the 
sum of the H-H distances on the horizontal axis, and their difference along the vertical axis. 
When no constraints were imposed, the equilibrium configuration was such that rf=0.81 A 
and z=1.97 A. As the total length of the molecule was shortened, both H-H distances 
gradually approached each other until equalization was observed when {d + z) fell below 
2.1 A. 

For comparison we also provide data extracted from lp-RbH 5 when it was reoptimized 
at intermittent pressures between 2-100 GPa. Between 17-20 GPa, equalization is observed 
when the total length of the H3" motifs nears 2.0 A. The onset of equalization occurs at very 
similar {z + d) values in both the molecule and in the solid. 

The lp-RbHg structure also has asymmetric H~---H 2 arrangements. However, equalization 
is not achieved at 100 GPa in spite of (d + z) falling below 2.0 A. We attribute this behavior 
to RbHg's structural peculiarities: whereas in RbH 5 each H~ has only one H 2 contact, in 
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FIG. 3: The difference between the z and d distances in an H~ • • • H2 fragment (see inset) is plotted 
vs. their sum for an isolated molecule in the gas phase. During the optimization, the total 
length of the molecule was constrained and the central atom was allowed to relax freely. These 
values are also illustrated for an H~ • • • H2 unit in lp-RbHg and lp-RbHs. These phases were 
reoptimized in order to determine the effect of pressure on (d + z) and {d — z). 

RbHg there are four. As a result, the H 2 unit is substantially shorter in RbH 9 than in RbH 5 
(it measures ^0.85 A in lp-RbH 5 , but never surpasses 0.78 A in RbH 9 even at 100 GPa). 
The total length of the fragment is substantially shorter at 2 GPa for RbH 5 than it is for 
RbH 9 suggesting that in the former phase the hydrogen sublattice is more compressed. 

If sufficiently squeezed the core electrons of the alkali metals start to overlap, the bands 
broaden and the valence electrons become localized in the interstitial regions; Na, for ex- 
ample, becomes an insulator at ^200 GPa [37J. For K, Rb and Cs, pressure is also known 
to induce an electronic s d transition, and by ^50 GPa Rb is in essence a d 1 metal 
[38]. Compress further, and the 5d band hybridizes with the broadened Ap semicore states. 
In Ref. [38J it was shown that both the core repulsion (a result of the overlap), and the 
hybridization are important factors in determining the most stable elemental Rb structures. 
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FIG. 4: (a) The bands of Cmcm hp-RbH 5 at 100 GPa calculated with a full Muffin Tin Orbitals of 
order N (iVMTO) basis are given in black. The green bands were computed with a minimal basis 
describing the Rb 4p, H2-0-, and HOMO-1 bands, and the red bands the HOMO. Along 
these high-symmetry lines the red and green bands do not cross. We illustrate the centered 
WFs obtained from the basis sets which yielded the (b) green, and (c) red bands. 

As illustrated in the Supporting Information (SI), the bands near the Fermi level in hp- 
RbH 5 at 100 GPa display some Rb 5d character. This is a minor contribution however, since 
the valence electron has been ionized into the hydrogen sublattice. Much more important 
for the electronic structure at this pressure is the overlap of the 4p semicore orbitals. In 
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fact, at 100 GPa even the Rb As orbitals interact with each other, as evidenced by their 
2.2 eV bandwidth. The substantial core overlap affects the electronic structure of hp-RbH 5 : 
a calculation of the DOS of the (H 5 )~ sublattice shows this hypothetical system to be a 
good metal (see the SI). 

Because of their overlap and hydridization with the Rb Ap bands, it is not possible to 
construct a Wannier function (WF) which describes only the H 2 -a bands. One can, however, 
construct a minimal basis set [39, HQ] that exactly reproduces the ten bands whose energies 
fall in between -18 to -5 eV in Fig. [i](a). The WFs obtained from this basis are the Rb p 
non-bonding and H 2 -a bonding (both not shown) , as well as the delocalized H3 a-bonding 
orbital illustrated in Fig. [3](b). One may also obtain a WF which yields just the two red 
valence bands. This WF, provided in Fig. [4](c), is primarily H3 non-bonding, but it also 
contains some Rb p and H 2 -a character. The filling of the two orbitals illustrated in Fig. 
[I](b, c) is a hallmark of a 3c-4e bond as would be found in the symmetric H3 , confirming 
our analysis of the bonding in hp-RbH 5 . 

As interesting as RbH 5 is, unfortunately it cannot be superconducting at least up to 
100 GPa where the PBE band gap is ~1.7 eV. In fact, the symmetrization of H3" is likely 
detrimental to the metallization. We calculate the HOMO-LUMO gap of an H3" molecule 
with z = d = 0.9Aas being 6.2 eV — over 40% larger than the 4.4 eV gap in the asymmetric 
fully optimized species. 

B. Metallization of RbH n 

The ^13 eV HOMO-LUMO gap in pure H 2 remains large as the molecules are brought 
together and bands form. The c^-bonding and <r*-antibonding bands undergo pressure- 
induced broadening, and are predicted to eventually overlap. But in experiments this does 
not happen below 342 GPa [T5] . 

One strategy to achieve metallization at experimentally achievable pressures is to perturb 
the electronic structure of H 2 by the addition of impurities [22] . Doping with an electroposi- 
tive element can lead to the partial filling of the H 2 a* bands and metallicity already at 1 atm 
as in LiH 6 [27]. Unfortunately, LiH 6 remains unstable with respect to decomposition into H 2 
and LiH below 110 GPa. Another possibility is the formation of an H~ donor-impurity band 
falling between the a u and a* bands as in LiH 2 . But the overlap of the Li Is core orbitals 
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hinders metallization, so this system remains a stable semi-metal above 120 GPa. Even 
though the sodium polyhydrides become viable at ^25 GPa, a metallic species which lies on 
the convex hull emerges only at 250 GPa. This NaH 9 phase consists of an Na + /H~ sublattice 
filled with H 2 molecular guests, and overlap of the Na 2p cores forestalls metallization. 

The pressure-dependance of the PBE band gaps of RbH 9 , RbH 5 , and RbH 3 (species 
which lie on the convex hull at some point within 2-250 GPa) are shown in Fig. [5j For all of 
these structures the calculated MIT transitions occur near 220 GPa, with the gap in RbH 9 
closing at slightly lower, and in RbH 5 at slightly higher pressures. 



i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 



> 

CD 

CL 
CD 

"D 
C 
CD 
CO 



3.0 
2.5 
2.0 
1.5 
1.0 
0.5 




■ 


RbHg 


• 


RbH 5 


A 


RbH3 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



25 50 75 100 125 150 175 200 225 250 

Pressure / GPa 



FIG. 5: The evolution of the PBE band gaps of the lowest-enthalpy RbHg, RbHs, and RbH3 
structures with increasing pressure. The most thermodynamically stable RbH n species, ie. the 
global minima (GM), are shown as unfilled symbols and are traced with a gray line. RbHg is also 
stable with respect to decomposition at 250 GPa, and is metallic throughout its domain of stability. 
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Two other structures, RbHy and RbH 8 , were found to metallize by 100 GPa. However, 
their enthalpies of formation fell above the convex hull and their calculated phonon 
dispersions showed soft modes. At 150 GPa P4/nmm-RbH 7 — with Rb + cations, H 2 
molecules, and H~ • • • H 2 moieties — was found to be mechanically stable, but it's AH F 
still lay above the convex hull. 

Structures of RbH 9 and RbH 5 near the MIT 

The behavior of compressed RbH 9 , which remains stable with respect to decomposition up 
to 250 GPa, is akin to NaHg. Both systems can be thought to consist of an M + /H~ lattice 
within a sea of H 2 molecules, and the band gap between the occupied hydridic bands and the 
H 2 a* bands eventually closes as the pressure is increased to ^200 GPa. P63/mmc-RbH 9 at 
150 GPa [56], see Fig. |6](a), was confirmed to be mechanically stable based on the absense 
of imaginary phonon modes. 

The total electronic densities of states (DOS) of this phase at 200 GPa is provided at the 
top of Fig. [6](d). The projected densities of states (see the SI) confirm the broadening of the 
Rb 4<s levels around -30 eV, the mixing of Rb 4p with the H 2 a u bands, and the contribution 
of the hydridic ions to the valence band. In addition, the H 2 -hydrogens furthest away from 
the H~ atoms contribute more to the valence DOS than their nearest neighbors. This type 
of density distribution is consistent with the HOMO in a system containing an assymmetric 
3c-4e bond (see the SI), confirming the interaction of the H 2 and H~ entities via delocalized 
multicentered bonding. 

At 220 GPa the enthalpy of Cm-RbH 9 recovered from an EA search at 250 GPa be- 
comes lower than that of the hexagonal phase (see Fig. [6](c)). Examination of the extended 
structure shown in Fig. [6](b) may lead one to classify it as an arrangement of Rb + cations, 
H 2 molecules, and slightly asymmetric H3 fragments. The second-nearest neighbor dis- 
tance about each hydrogen atom in Cm-RbH 9 lies between 0.93-1.16 A, as compared with 
1.27-1.38 A in the P63/mmc-symmetry structure at 220 GPa. The more tightly-packed 
Cm-phase has a slightly lower volume, by ^2% at 220 GPa, and is therefore preferred at 
higher compression. This phase is metallic with a much higher density of states at the Fermi 
level than the hexagonal structure. Such a scenario is reminiscent of the sodium polyhy- 
drides for which the enthalpy of a metallic Cmc2i-NaH 9 structure became lower than that 
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FIG. 6: Supercells of the lowest enthalpy RbHg phases recovered at (a) 150 GPa (P6s/mmc- 
symmetry) , and (b) 250 GPa (Cm-symmetry); and (c) their enthalpies of formation at various 
pressures, (d) Total densities of states of (top) P63 / 'mmc-RbHg at 200 GPa, and (bottom) Cm- 
RbHg at 220 GPa. The Fermi levels are set to zero. 

of insulating Pm -NaH 9 at 300 GPa. However, due to geometry convergence issues, we were 
unable to confirm the mechanical stability for Cm-RbH 9 . 

The Cmcm hp-RbH 5 phase illustrated in Fig. [2](a) and reproduced in Fig. [7](a) remains 
the most stable RbH 5 structure up to 250 GPa. In going from 100 to 250 GPa, its band gap 
decreases from to 1.7 to 0.2 eV, see Fig. [7](d). The convex hull, Fig. [jjc), shows that RbH 5 
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FIG. 7: Supercells of the lowest enthalpy RbHs phases recovered at (a) 100 GPa (Cmcm- 
symmetry) , and (b) 250 GPa (Cmmm-symmetry); and (c) their enthalpies of formation at various 
pressures, (d) Total densities of states of (top) Cmcm-KbB.^ at 100 GPa, (middle) Cmcm-KbB.^ 
at 250 GPa, and (bottom) Cmmm-RbHs at 250 GPa. The Fermi levels are set to zero. 



becomes slightly unstable with respect to decomposition into RbH 6 and RbH 3 at 250 GPa. 
The EA search at 250 GPa also recovered a metallic Cmmm-RbH 5 arrangement shown in 
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Fig. [fjb). The two structures possess very similar interatomic distances, and are within 
dynamical enthalpies of each other (2 meV/atom), in stark contrast to the situation at 
100 GPa where Cmcm-RbH 5 was at least 15 meV/atom lower in enthalpy than any other 
structure recovered in the EA searches. The zero point energy (ZPE) corrections favor 
the mechanically stable Cmmm structure by ^5 meV/atom over the Cmcm at 200 GPa, 
causing the ZPE-corrected enthalpies to fall within 5 meV/atom. The closeness of the 
enthalpies under pressure suggests the potential melting of the hydrogen sublattice, and 
quantum fluid-like behavior [41J. 



RbH 3 and RbH 6 : H3 molecules and their eventual polymerization 

RbH 3 , which becomes stable with respect to decomposition at 100 GPa, does not exhibit any 
imaginary phonon modes at 150 GPa [57] . By 220 GPa, this 74i/amd~symmetry structure 
(shown in Fig. [8](a)) overtakes RbH 5 to become the most stable polyhydride. Its reign, 
however, is shortlived as the (mechanically stable) Cmmm-RbH 3 structure in Fig. [8](b) 
becomes enthalpically preferred at 220 GPa [58], see Fig. [8](c). The difference between the 
ZPE corrections for the two phases was negligible at 220 GPa (within 1 meV/atom). 

Both RbH 3 configurations are composed of Rb + cations and H3 anions spread about a face 
centered cubic lattice, which is distorted somewhat in the IAi/amd structure. The shortest 
intra- and inter-Hg H-H distances are quite similar in the two, measuring 0.87 and 1.47 A, 
respectively, at 220 GPa. The biggest qualitative difference between them is the stacking 
within the lattice along the axis of the H3 units: an alternating/repeating sequence of Rb + 
and H3 is found in the 7*4i/ 'amd/Cmmm arrangements. The emergent layering of Rb + and 
H3 in the latter is shown in Fig. [8](b), wherein the layers lie in a plane perpendicular to the 
page. The Cmmm structure is metallic at 220 GPa, whereas the IAi/amd is a semi-metal 
as shown in Fig. [8](d). 

At 250 GPa, both RbH 3 and RbHg are still stable but the convex hull for RbH n with 
n=3-9 flattens to the extent that Imma-RbH 6 shown in Fig. [9](a) becomes viable as well 
[59] . Interestingly, the hydrogen sublattice of the RbH 6 phase can no longer be thought of 
as being made up of molecular fragments, but instead of arrangements of H3 units brought 
together to form what are essentially chains, ie. (Hg - )^. At this pressure RbH 6 , whose DOS is 
illustrated in Fig. Mb)^ is a good metal. Whereas the HOMO of a symmetric H3 is localized 
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Energy / eV 

FIG. 8: Supercells of the lowest enthalpy RbH3 phases recovered at (a) 150 GPa {l\\jamd- 
symmetry) , and (b) 250 GPa (Cramm-symmetry); and (c) their enthalpies of formation at various 
pressures, (d) Total densities of states of (top) /4i/amd-RbH3 at 220 GPa, and (bottom) Cmmm- 
RbHs at 220 GPa. The Fermi levels are set to zero. 

on the two ends of the molecule (see the SI), every hydrogen atom contributes about the 
same amount to the projected DOS just below the Fermi level in RbH 6 , suggesting that a 
more delocalized valence electron distribution is forming along (H3 ) oc . 

Phonon calculations (in the harmonic approximation) on RbH 6 suggest the structure is 
unstable — a couple of large-amplitude imaginary frequencies are found even at the T-point. 
Visual inspection of the soft modes indicates that they correspond to collective motions 
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FIG. 9: (a) Supercell, and (b) total density of states of the lowest enthalpy Imma-RbHQ phase 
recovered at 250 GPa. The Fermi level is set to zero. 

of the hydrogen atoms along the chains. To investigate this further, molecular dynamics 
(MD) simulations were run on a supercell of Imma-RbH 6 . The evolution of a number of 
interatomic H-H distances from the 250 GPa-optimized RbH 6 structure throughout a 5 ps 
NVT MD simulation is shown in Fig.[l0|(a). The trajectories are shaded to illustrate whether 
the H-H distances corresponded to those along the constituent "chains" (black) , or between 
them (gray). A similar plot is provided for RbH 3 , where the intramolecular H-H distances 
are all slightly under 1.0 A, and the intermolecular ones measure 1.5 A or higher. These 



distances oscillate only slightly about their equilibrium value, see Fig. 10(b). In RbH 6 , on 



the other hand, the black and gray lines overlap at about 1.1 A, suggesting that the hydrogen 
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sublattice is polymerizing. 
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FIG. 10: The evolution of the H-H distances within the (a) (Hg ) OQ "chains" and the (b) Hg~ 
molecules (black lines) , and between them (gray lines) during a 5 ps molecular dynamics simulation 
of (a) Imma-RbHg and (b) Cmmm-RhYi^ at 250 GPa. 

Optimized snapshots from the MD simulation, wherein a "bond" has been drawn between 
two hydrogen atoms if their interatomic distance measures less than 1.2 A, are provided in 
Fig. [TT|a). At first sight it may appear that each configuration maintained the distinct 
chain-like motifs in the original structure, but the main difference between them lies in the 
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FIG. 11: (a) Optimized structures of snapshots taken from a 5 ps molecular dynamics simulation 
of Imma-RbHg at 250 GPa, and (b) the initial and final structures. In (b), six hydrogen atoms, 
and their images, are enlarged and highlighted in a 1 x 2 x 2 represent ion of the simulation cell to 
show their displacement over the course of the simulation. 
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relative positions of the hydrogen atoms within the chains. For example, at t = 5 fs they are 
eclipsed along the axis perpendicular to the page, whereas at 3 and 5 ps they are staggered. 
All three of the optimized structures shown in Fig. [TT](a) are within dynamical enthalpies 
(3 meV/atom) of each other. 

The presence of multiple energetically degenerate local minima and the thermally 
accessible barriers between them (as indicated by the short molecular dynamics run) 
suggest the hydrogen sublattice may even exhibit liquid-like behavior [41], akin to that 
predicted for the high-pressure forms of PbH 4 [42j . This is further corroborated by noting 
the rise in inter-chain H-H distances toward the end of the MD-simulation in Fig. [lO|(a) is 
caused not by chains moving far apart from each other, but by the increase in interatomic 
distance between two atoms on adjacent chains as they flow. In Fig. [TT|(b), the initial and 
final positions of select hydrogen atoms are enlarged to show their movement during the 
simulation. On average, each hydrogen atom moved 1.9 A (with a standard deviation, 
a, of 1.2) while each Rb atom moved only 0.07 A (a=0.02). Furthermore, the distances 
between and within the H3 building blocks of the (H3 ) 00 chains (the solid vs. dashed lines 



in Fig. 10(a)) are indistinguishable from one another. Since these distances measure 0.99 A 
and 0.87 A in the original Imma-RbH 6 structure, only a 12% difference, this suggests 
that the thermal vibrations of the H atoms about the chain are enough to circumvent 
the distinction of the two characteristic H-H distances. The hydrogenic sublattice in this 
structure can be thought to consist of flowing, partially negatively charged, polymeric chains. 



III. CONCLUSIONS 

Our first-principles computations predict a number of rubidium poly hydrides, RbH n 
n > 1 , which become stable at pressures as low as 2 GPa. There is a wide variety in their 
structures, and their hydrogenic sublattices can be viewed as being composed of H~, H 2 and 
H3 units. We note the similarities with polyiodide compounds whose complex, and often 
times polymeric-like lattices can be understood in terms of interactions between the I - , I 2 , 
as well as the symmetric and assymmetric linear (or nearly) I3 building blocks [43] . 

RbH 5 , the most stable structure between 15-220 GPa contains Rb + and H 2 molecules, 
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along with the symmetric H3 motifs — the simplest example of a 3c-4e bond. The broadened 
Rb Ap bands hybridize with hydrogenic bands, and the overlap of the Rb As/Ap orbit als makes 
this phase an insulator. Compress further and a metallic RbH 3 phase made up Rb + and H3 
building blocks becomes the most stable poly hydride. At 250 GPa, alternative forms of the 
hydrogen sublattice emerge in RbH 6 in the form of polymeric chains. Molecular dynamics 
simulations suggest the hydrogen sublattice in RbH 6 may exhibit one dimensional liquid-like 
behavior. 



IV. COMPUTATIONAL DETAILS 



The structural searches were performed using the open-source EA XtalOpt Release 7 



(http://xtalopt.openmolecules.net/), and the parameter set suggested in Ref. [29]. EA 
searches were run on simulation cells with two RbH n formula units at 2, 10, 30, 100, 150 
and 250 GPa. The lowest enthalpy structures from each search were relaxed in the pressure 
range from 0-250 GPa. 

First principles DFT calculations were performed with the Vienna ab initio Simulation 
Package (VASP) version 4.6.31 [44]. The projector augmented wave (PAW) method [45j was 
used to treat the core states, and a plane-wave basis set with an energy cutoff of 600 eV 
was employed. An energy cutoff of 500 eV was used in the structural searches. The Rb 
4<s/4p/5<s electrons were treated explicitly in all of the calculations. The gradient-corrected 
exchange and correlation functional of Perdew-Burke-Ernzerhof (PBE) [46] was adopted, 
and the fc-point grids were generated using the T-centered Monkhorst-Pack scheme. The 
number of divisions along each reciprocal lattice vector was chosen such that the product 
of this number with the real lattice constant was 20 A in the structural searches and 40 A 
otherwise. 

Phonons and thermodynamic properties of select RbH 3 , RbH 5 , RbH 6 , and RbH 9 phases 
were calculated using the PHON [47] package. The supercells used for the phonon calcula- 
tions were chosen such that the number atoms in the simulation cell was always between 96 
and 216 atoms. We have neglected the zero point energy (ZPE) corrections to the enthalpy 
in most of our discussion. Since the importance of ZPE contributions when taking into ac- 
count the stability of systems with light atoms has been well-documented in the literature 
[32] . we provide two examples in the text of how the inclusion of the ZPE term influences the 
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stability RbH 3 and RbH 5 at 200-220 GPa. Overall, our assessment is that at the pressures 
we focus on, the differences in ZPE between RbH n structures are never large. 

In order to assess the mechanical stability of RbH 6 , ab-initio molecular dynamics (MD) 
simulations were run on RbH 6 using VASP. The MD simulations were performed under 
constant volume and temperature conditions (NVT) for 5 ps with 0.5 fs timesteps. The 
initial velocities were initialized randomly according to a Maxwell-Boltzmann distribution 
such that the starting temperature of the system was 400 K and the temperature was 
maintained through intermittent velocity rescaling. A 96-atom simulation cell was used for 
RbH 6 and a 128-atom simulation cell for the complementary RbH 3 simulation. 

The band structures of select phases were calculated using the tight-binding linear muffin- 
tin orbital (TB-LMTO) method @Hl |49], and for hp-RbH 5 NMTO methods (MTOs of 
order N) [391 150] . In the TB-LMTO calculations the VWN [51] local exchange correlation 
potential was used along with the Per dew- Wang [52] generalized gradient approximation. 
Scalar relativistic effects were included. The structural parameters of hp-RbH 5 were taken 
from those optimized with VASP at 100 GPa, and four different types of empty interstitial 
spheres were inserted, each with a single s function. The calculations utilized 2768 irreducible 
points in the tetrahedron k— space integrations [53]. The TB-LMTO densities of states 
showed good agreement to those obtained with VASP (see the SI). The fat bands, whose 
width is proportional to a particular orbital character, are also provided in the SI. 

The band structure with a full AHV1TO basis (Rb 4s, 4p, 5d and H Is, 2p) of RbH 5 at 
100 GPa corresponded well with the bands of Fig. [I](a). The a Rb 4p + H 2 a + H3 bonding" 
bands in this figure were computed with an A^MTO minimal basis consisting of s and p 
functions on all Rb atoms, an s function on every second H 2 atom, and an s function on 
the central H atom in an H3 unit. The "Hg non^bonding" bands were computed with an s 
function on one of the terminal hydrogens in each H3 unit. In the SI, we show these bands 
along with the energy meshes used in the calculations. 



22 



Acknowledgements 



We acknowledge the NSF (DMR- 1005413) for financial support, and the Center for Com- 
putational Research at SUNY Buffalo for computational support. 



[1] J. J. Thomson, Phil. Mag., 1912, 24, 209-253. 

[2] R. Golser, H. Gnaser, W. Kutschera, A. Priller, P. Steier, A. Wallner, M. Cizek, J. Horacek 

and W. Domcke, Phys. Rev. Lett, 2005, 94, 223003 (1-4). 
[3] H. Gnaser and R. Golser, Phys. Rev. A, 2006, 73, 021202 (1-4). 

[4] W. Wang, A. K. Belyaev, Y. Xu, A. Zhu, C. Xiao and X. F. Yang, Chem. Phys. Lett., 2003, 
377, 512-518. 

[5] C. A. Coulson, Proc. Camb. Phil. Soc, 1935, 31, 244. 
[6] M. E. Schwartz and L. J. Schaad, J. Chem. Phys., 1967, 47, 5325-5334. 
[7] D. Stevenson and J. Hirschfelder, J. Chem. Phys., 1937, 5, 933-940. 
[8] J. Starck and W. Meyer, Chem. Phys., 1993, 176, 83-95. 
[9] F. Robicheaux, Phys. Rev. A, 1999, 60, 1706-1709. 
[10] M. Ayouz, O. Dulieu, R. Guerout, J. Robert and V. Kokoouline, J. Chem. Phys., 2010, 132, 
194309 (1-11). 

[11] A. K. Belyaev, A. S. Tiukanov and W. Domcke, Chem. Phys., 2006, 325, 378-388. 

[12] N. W. Ashcroft, Phys. Rev. Lett, 1968, 21, 1748-1749. 

[13] N. W. Ashcroft, Phys. Rev. Lett, 2004, 92, 187002 (1-4). 

[14] P. Loubeyre, F. Occelli and R. LeToullec, Nature, 2002, 416, 613-617. 

[15] C. Narayana, H. Luo, J. Orloff and A. L. Ruoff, Nature, 1998, 393, 46-49. 

[16] P. Cudazzo, G. Profeta, A. Sanna, A. Floris, A. Continenza, S. Massidda and E. K. U. Gross, 

Phys. Rev. Lett, 2008, 100, 257001 (1-4). 
[17] G. Gao, A. R. Oganov, A. Bergara, M. Martinez-Canales, T. Cui, T. Iitaka, Y. Ma and G. Zou, 

Phys. Rev. Lett, 2008, 101, 107002 (1-4). 
[18] J. S. Tse, Y. Yao and K. Tanaka, Phys. Rev. Lett, 2007, 98, 117004 (1-4). 
[19] C. J. Pickard and R. J. Needs, Phys. Rev. B, 2007, 76, 144114 (1-5). 

[20] I. Goncharenko, M. I. Eremets, M. Hanfland, J. S. Tse, M. Amboage, Y. Yao and I. A. Trojan, 



23 



Phys. Rev. Lett, 2008, 100, 045504 (1-4). 
[21] G. Gao, H. Wang, A. Bergara, Y. Li, G. Liu and Y. Ma, Phys. Rev. B, 2011, 84, 064118. 
[22] D. D. Klug and Y. Yao, Phys. Chem. Chem. Phys., 2011, ASAP, doi: 10.1039/clcp21791g. 
[23] J. Feng, R. G. Hennig, N. W. Ashcroft and R. Hoffmann, Nature, 2008, 451, 445-448. 
[24] M. Somayazulu, P. Dera, A. F. Goncharov, S. A. Gramsch, P. Liermann, W. Yang, Z. Liu, 

H. Mao and R. J. Hemley, Nat. Chem., 2010, 2, 50-53. 
[25] T. A. Strobel, M. Somayazulu and R. J. Hemley, Phys. Rev. Lett, 2009, 103, 065701 (1-4). 
[26] S. Wang, H. Mao, X. J. Chen and W. L. Mao, Proc. Nat. Acad. Set, 2009, 106, 14763-14767. 
[27] E. Zurek, R. Hoffmann, N. W. Ashcroft, A. R. Oganov and A. O. Lyakhov, Proc. Natl. Acad. 

Sci., 2009, 42, 17640-17643. 
[28] P. Baettig and E. Zurek, Phys. Rev. Lett, 2011, 106, 237002 (1-4). 
[29] D. C. Lonie and E. Zurek, Comput. Phys. Commun., 2011, 182, 372-387. 
[30] D. Lonie and E. Zurek, Comput. Phys. Commun., 2011, 182, 2305-2306. 
[31] R. Ahuja, O. Eriksson and B. Johansson, Physica B, 1999, 265, 87-91. 
[32] C. J. Pickard and R. J. Needs, Nature Phys., 2007, 3, 473-476. 

[33] W. Grochala, R. Hoffmann, J. Feng and N. W. Ashcroft, Angew. Chem. Int. Ed., 2007, 46, 
3620-3642. 

[34] J. M. McMahon and D. M. Ceperley, Phys. Rev. Lett, 2011, 106, 165302 (1-4). 

[35] V.-F. Tebbe and U. Georgy, Acta. Cryst, 1986, C42, 1675. 

[36] F. Maseras, A. Lledos, E. Clot and O. Eisenstein, Chem. Rev., 2000, 100, 601. 

[37] Y. M. Ma, M. Eremets, A. R. Oganov, Y. Xie, I. Trojan, S. Medvedev, A. O. Lyakhov, M. Valle 

and V. Prakapenka, Nature, 2009, 458, 182-185. 
[38] A. McMahan, Phys. Rev. B, 1984, 29, 5982. 

[39] E. Zurek, O. Jepsen and O. K. Andersen, ChemPhysChem, 2005, 6, 1934-1942. 

[40] E. Zurek, O. Jepsen and O. K. Andersen, Inorg. Chem., 2010, 49, 1384-1396. 

[41] S. A. Bonev, E. Schwegler, T. Ogitsu and G. Galli, Nature, 2004, 431, 669-672. 

[42] P. Zaleski-Ejgierd, R. Hoffmann and N. W. Ashcroft, Phys. Rev. Lett, 2011, 107, 037002. 

[43] P. H. Svensson and L. Kloo, Chem. Rev., 2003, 103, 1649-1684. 

[44] G. Kresse and J. Hafner, Phys. Rev. B, 1993, 47, 558-561. 

[45] P. Blochl, Phys. Rev. B, 1994, 50, 17953. 

[46] J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett, 1996, 77, 3865-3868. 

24 



[47] D. Alfe, G. D. Price and M. J. Gillan, Phys. Rev. B, 2001, 65, 045123. 

[48] O. K. Andersen, Phys. Rev. B., 1975, 12, 3060. 

[49] O. K. Andersen and O. Jepsen, Phys. Rev. Lett, 1984, 53, 2571. 

[50] O. K. Andersen and T. Saha-Dasgupta, Phys. Rev. B., 2000, 62, R16219. 

[51] S. H. Vosko, L. Wilk and M. Nusair, Can. J. Phys., 1980, 58, 1200. 

[52] J. P. Perdew and Y. Wang, Phys. Rev. B., 1986, 33, 8800. 

[53] P. E. Blochl, O. Jepsen and O. K. Andersen, Phys. Rev. B., 1994, 49, 16223. 

[54] At 100 GPa: a=2.97 A, 6=7.39 A, c=4.88 A with the atoms at the following Wyckoff sites: H 

4/ (0.500 0.781 0.171), H 4/ (0.500 0.360 0.064), H 2c (0.500 0.633 0.750), and Rb 2c (0.500 

0.923 0.750) 

[55] The molecular calculations on were performed using the ADF (www.scm.com) package 

(BP functional with DZP basis set ) 
[56] At 150 GPa: a=3.65 A, c=5.42 A with the atoms at the following Wyckoff sites: H 12k (0.535 

0.465 0.057), H 4/ (0.333 0.667 0.317), H 2d (0.667 0.333 0.250), and Rb 2b (0.000 0.000 0.250) 
[57] At 150 GPa: a=2.83 A, c=9.70 A with the atoms at the following Wyckoff sites: H 4e (0.500 

0.250 0.217), H 2b (0.500 0.250 0.125), and Rb 2a (0.000 0.250 0.875) 
[58] At 220 GPa: a=4.65 A, 6=2.71 A, c=2.70 A with the atoms at the following Wyckoff sites: 

H 2g (0.188 0.000 0.000), H la (0.000 0.000 0.000), and Rb lc (0.500 0.000 0.500) 
[59] At 250 GPa: a=6.54 A, 6=3.70 A, c=4.75 A with the atoms at the following Wyckoff sites: H 

4/ (0.244 0.747 0.155), H 4/ (0.099 0.391 0.485), H 4/ (0.888 0.607 0.794), and Rb 2e (0.000 

0.121 0.750) 



25 



Table of Contents Text: 

A pressing motif: At 15 GPa RbH 5 is predicted to become stable with respect to de- 
composition into H 2 and RbH. The unit cell contains Rb + cations, H 2 molecules, and the 
symmetric linear H3 unit — the simplest example of a three center four-electron bond. At 
higher pressures H3 is further stabilized in an RbH 3 phase which metallizes at 220 GPa. 
Keywords: hydrides, alkali metals, ab-initio calculations, high-pressure chemistry, bond 
theory 
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